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We discuss the structural and electronic properties of tetragonal CuO grown on SrTiO3(100) 
by means of hybrid density functional theory. Our analysis explains the anomalously large Cu-0 
vertical distance observed in the experiments (?^2.7A) in terms of a peculiar frustration between 
two competing local Cu-0 environments characterized by different in-plane and out-of-plane bond 
lengths and Cu electronic populations. The proper inclusion of substrate effects is crucial to under- 
stand the tetragonal expansion and to reproduce correctly the measured valence band spectrum for 
a CuO thickness of 3-3.5 unit cells, in agreement with the experimentally estimated thickness. 

PACS numbers: 68.35.bt, 73.61. Le, 71.15.-m 
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I. INTRODUCTION 

Very recently, Siemons et alX synthesized CuO with 
a tetragonal (i.e. elongated rock salt structure) on 
a SrTiO3(100) substrate. In a subsequent theoretical 
study^ we investigated bulk phases of tetragonal CuO 
by applying a hybrid density functional theory approach, 
which is able to deal with the electronic, ionic and mag- 
netic properties of such a system. Studying a variety of 
magnetic orderings we obtained two energy minima as a 
function of tetragonal distortion. The antiferromagneti- 
cally ordered phase (denoted as TET2) with a tetragonal 
distortion of c/a=1.377 appeared to be the most stable 
one with its lattice parameters being in very good agree- 
ment with the experimental results. Furthermore, for the 
TET2 phase we predicted a very high Neel temperature 
(Tat) of 800 K, which would perfectly fit the trend of Tat 
for the 3d-transition metal monoxides^. Also, the calcu- 
lated density of states of the TET2 phase agreed well with 
the experimental valence band spectrum with the excep- 
tion of a few residual differences related to the structure 
of the main peak and the presence of an additional peak 
in the low energy spectrum which we have attributed to 
surface/substrate effects not taken into account in our 
previous analysis. 

The aim of the present work is now, to corroborate 
our basic findings for the artificial TET2 phase by mod- 
elling experiment as realistically as possible. This is what 
we do by describing the adsorbate system in terms of 
an (CuO)[n]/SrTiO3(100) slab for varying thicknesses of 
CuO layers [with n denoting a CuO coverage of 2,3,3.5 
and 4 unit cells (u.c.)]. 

Shortly after the publication of our bulk CuO study, 
a further theoretical study appeared, which was based 
on the self-interaction-corrected local-density functional 
method^ Again, for a variety of magnetic orderings two 
energy minima were found as a function of the tetrag- 
onal distortion, but the c/a values were much closer to 



1 than in our work and, consequently, in worse agree- 
ment with the experimentally observed c/a ratio of 1.357. 
Therefore, the present study also serves the purpose to 
buttress the theoretical capability of our applied hybrid 
density functional theory approach and to understand 
the nature of the exceptionally large vertical elongation 
observed in SrTiOa supported CuO. 



II. COMPUTATIONAL ASPECTS 

Because standard density functional theory applica- 
tions fail in correctly describing the ground state of 3d- 
transition metal monoxides, a more sophisticated (and 
much more costly) approach has to be chosen. Therefore, 
for the present study we apply hybrid density functional 
theor}^ based on the Heyd-Scuseria-Ernzerhof (HSE) 
method^ as implemented in the Vienna ah initio sim- 
ulation package (VASP)^— . All the corresponding tech- 
nical parameters were the same as in our recent study^ of 
tetragonal bulk CuO. In particular, 1/4 of short-ranged 
Hartree Fock exchange was admixed to the generalized- 
gradient-approximation exchange-correlation functional. 
The adsorbate systems (CuO)[n]/SrTiO3(100) (n=2,3,3.5 
and 4 unit cells) were modelled by a repeated slab scheme 
containing up to 62 atoms per slab for the highest cover- 
age. Figures [1] and [2] sketches the basic layer- wise atomic 
arrangements. All atomic positions of CuO layers were 
fully relaxed, whereas for the SrTiOs substrate the lowest 
5 bottom layers, corresponding to a full SrTiOa unit cell, 
were kept fixed. The lateral lattice parameter a=3.90 A 
for SrTiOs was taken from the recent VASP-based HSE 
study by R. Wahl and coworkers^. Indeed, the optimal 
HSE value of a is in very good agreement with the mea- 
sured lattice constant aExpt=3.900 A ^^. In order to make 
the HSE computations feasible, ferromagnetic (FM) or- 
dering was assumed although a specific antiferromagnetic 
ordering appeared to be energetically more stable in the 
ideal tetragonal bulk phase. ^ Nevertheless, also because 
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FIG. 1: (Color online) Top view of the interface models of 
the (CuO)[n]/SrTiO3(100) adsorbate system as studied by the 
present HSE approach. The numbers in the squares (bottom 
left corner of each structure) denote the calculated energy 
difference (in meV/A^) relative to the Ti02 terminated Ml 
stacking, which is the most stable arrangement. 



of the small magnetic Cu- mo merits of about 0.7 jj^b we 
expect only a small influence of the magnetic ordering on 
the structural relaxations and valence band spectrum^-, 
which are the aim of study of our present work on the 
(CuO)[n]/SrTiO3(100) substrate system. The k-points 
integration for the structural optimization runs has been 
done by using a 4x4 two dimensional Monkhorst-Pack 
grid, which was increased to 6x6 for the final electronic 
relaxation. 



III. RESULTS AND DISCUSSION 
A. Structural Properties 

For finding the energetically most stable atomic ar- 
rangement several terminations of the SrTiO3(100) sub- 
strate and stacking of the CuO adsorbate have been in- 
vestigated. Fig. [1] sketches all the studied cases: the 
substrate is either Ti02 or SrO terminated. CuO layers 
may be accommodated in several ways: (i) Cu on top 
of oxygen atoms (Ml), (ii) Cu on top of Ti or Sr atoms 
(M2) or (iii) Cu and O in hollow sites (M3). These (FM) 
calculations have been done by placing 2 CuO unit cells 
(5 CuO layers) on SrTiOs (6 unit cells). Fig. [1] shows 
that the Ti02 terminated Ml structure is the most stable 
one. It should, however, be noted that rather close in en- 
ergy (only 14 meV/A^ less stable) is the SrO terminated 
M2 stacking. 

From now on, we only discuss the Ti02 terminated Ml 
structure. We study the evolution of the structural and 
electronic properties of the (CuO)[n]/SrTiO3(100) adsor- 
bate system for a CuO coverage of n=2,3,3.5 and 4 unit 
cells. This choice is based on the experiments estimate, 
that a tetragonal CuO-layer consisting of 3 to 4 unit 
cells can be grown on SrTiOs (100)i . The corresponding 
structural models are displayed in Fig [2l 

The structural properties are summarized in Fig l3] and 
TablelH Fig. [3] illustrates that theoretical and experimen- 
tal geometrical properties agree well when in the calcu- 
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FIG. 2: (Color online) Side view of the (CuO)[n]/SrTiO3(100) 
slabs for a coverage by n=2,3,4 unit cells of CuO as used in 
the HSE calculation. (The n=3.5 coverage, which was also 
studied, is not shown). Atomic arrangements according to 
the Ti02 terminated substrate and the Ml stacking of the 
CuO adsorbate (see Fig. [T|). The layers are numbered with 
respect to the interface: I+j denotes CuO layers, I — j marks 
SrTiOs layers. The side projected two dimensional unit cell 
with lattice parameters a and c is sketched by the dashed 
rectangle. Detailed structural data are listed in Table H] 



lation the CuO coverage is 3.5 unit cells, which is also in 
the estimated range of the experimentally derived thick- 
ness (parameter D in Table H]). The experimental value 
of c/a= 1.357 is close to the theoretical value of 1.345, 
which compares favorably to c/a= 1.377 as obtained from 
the HSE study for the ideal bulk TET2 phase.^ In gen- 
eral, when increasing the number of unit cells c and c/a 
decrease as shown in Table [H 

The analysis of the interlayer distances reveals that in- 
side the SrTiOs substrate the interface effects are healed 
out rather rapidly. Only the distance of the layer closest 
to the interface (<ii-i/i-2) experiences a distinct short- 
ening, as it is also the case for the CuO layer at the 
interface (<ii+i/i-i). On the other hand, the CuO layer 
distances appear to be quite sensitive to the thickness 
of adsorbed CuO. For coverage of 2 and 3 unit cells the 
distances are significantly larger than the reference value 
of 2.69 A which was derived from the ideal TET2 bulk 
study. For coverage larger than 3 unit cells the distances 
between the inner CuO layers are significantly shortened, 
whereas the top 2 layers expand outwards. This change 
in bond length is reflected by the change of ionicity of 
Cu, as will be discussed in the following. 



B. Electronic Properties 

Figure m compares the HSE calculated density of states 
(DOS) with the experimental valence band spectrum. 
It is obvious that the total as well as the local CuO- 
projected DOS is rather sensitive to the thickness of the 
CuO block. Strikingly, the calculated DOS for the ideal 
TET2 bulk phase and for the system with n=3-3.5 unit 



TABLE I: Geometrical data for (CuO)[n]/SrTiO3(100) for 
n=2, 3, 3.5 and 4 (n: number of CuO unit cells). The data for 
c and c/a are based on averages of the vertical lattice param- 
eter according to Ref. 1, and the lateral parameter a=3.90 
A of SrTiO3(100) is taken from another VASP-based HSE 
stud}^. The total thickness of adsorbed CuO is described by 
D. The inter layer distance between two subsequent layers as 
illustrated in Fig. [2] is denoted by d, and (ii+i/i_i indicates 
the interface separation between SrTiOs and CuO. Layers 
are labelled according to FigO The lattice parameters of the 
HSE calculation for the ideal TET2 phase^ are a=3.908 and 
c=5.381 A resulting in an ideal bulk-like layer distance of 
2.69 A. Experimental values are included when available. All 
length parameters are given in units of A. 
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cells coverage agree well with experiment, quite in con- 
trast to all cases. Though the superstructure with 2 CuO 
u.c. already displays some film features which are not re- 
produced in the calculated bulk spectrum (such as the 
appearance of the peak A), the overall comparison with 
the measured curve is not satisfactory due to (i) a much 
too dominant presence of substrate-related states which 
weaken the intensity of the main CuO peaks around -4 eV 
and (ii) the lack of a one-to-one correspondence between 
theoretical and experimental peaks. 

As was discussed for the bulk phase TET2^, the agree- 
ment with experiment is rather good with the exception 
of one missing peak (indicated by the letter A in Fig. 2]) 
and one extra peak (indicated by the cross) in the calcu- 
lated spectrum. From the DOS for n=3-3.5 one clearly 
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FIG. 3: (Color online) Thickness dependent variation of the 
structural parameters D, c/a and c of the CuO adsorbate on 
SrTiO3(100) as listed in Table [D Horizontal bars illustrate 
the experimental data of Ref[l|. 



deduces that the experimental structure is reproduced 
very well (especially for 3.5 u.c.) and the overall agree- 
ment is improved with respect to the ideal bulk phase. 
In particular, a new peak (A) is found due to the CuO 
interface layer (the local CuO DOS has a very distinctive 
peak at this energy position, which is also reflected in 
the total DOS) and the spurious peak found in the ideal 
bulk phase disappears. Finally, according to the DOS for 
n=3.5 in Fig. |4]the highest experimental peak, can now 
be mainly attributed to the substrate, because there the 
local CuO DOS shows a weak depression, again in line 
with the observed valence band spectrum. 

Figure [5] illustrates the layer dependent change of ion- 
icity of the Cu atoms decomposed over the orbital quan- 
tum number I. A comparison with the corresponding 
bulk values for both the monoclinic and tenorite struc- 
tures is also given. Leaving out the layer distance at 
the interface there is a uniform trend to be seen: with 
increasing distance from the interface the Cu charge is 
reduced by about 0.15 e~ (i.e. the ionicity of the pos- 
itive Cu-ions increases). This change in charge/valency 
(topmost panel) is mostly due to the change of p- and 
d-like charges (2nd and 3rd panel from the top). For the 
two Cu layers closest to the surface the total Cu charge 
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FIG. 4: (Color online) Results of the HSE calculations: total 
(light grey) and local CuO-projected (dark green) density of 
states for (CuO)[n]/SrTi03 [n=2, 3, 3.5 and 4 unit cells (u.c.)] 
and for the ideal bulk tetragonal TET2 phase.— The experi- 
mental curve is the valence-band spectrum of Ref[l|. The cir- 
cles and crosses indicate the position of the measured peaks. 
Energy zero is chosen to be the top of the valence band. Posi- 
tion of the maximum of the experimental spectrum is adjusted 
to the corresponding maximum of the HSE density of states. 
The calculated DOS is broadened by a Gaussian of half- width 
0.2 eV. 



is even less than for the tetragonal bulk phase. This in- 
crease in valency is directly accompanied by an increase 
in the inter layer distance, as clearly illustrated by Table 

HI 

In general, an increase in bond length corresponds to 
an increase in ionicity: typical Cu-0 bond lengths involv- 
ing Cu^+ ions are of the order of 2 A whereas for Cu^+ 
the value might strongly increase up to 2.75 A (see Ref. 
[12I ). In our case, the increase in interlayer distance corre- 
sponds to the stretching of the Cu-0 bond in the vertical 
direction (perpendicular to the Cu-0 layers). Accord- 
ing to Table |I] for the topmost layer of the n=3.5 case 
the bond length is 2.87 A, being well in the bond length 
regime of Cu^+ ionicity. On the other hand, the in-plane 
bond lengths -due to the SrTiOs substrate or the lattice 
parameter a=3.91 A of the optimized bulk TET2 phase- 



correspond rather well to an ionicity of Cu^+. On the ba- 
sis of our present surface and recent bulk TET2 study it 
seems that the Cu atoms appear in two competing con- 
current ionic states, depending on the directions of the 
Cu-0 bonds. 

Our interpretation of the layer-dependent modulation 
of the Cu electronic population and Cu-0 bond lengths 
suggesting a peculiar coexistence of Cu^+-like and Cu^+- 
like behavior is in line with the combined DFT and pho- 
toelectron spectroscopy study on copper oxide clusters 
presented by L.-S. Wang and coworker^. By investi- 
gating the structural and electronic properties of Cu20aj 
(x=l-4) these authors report a significant change both in 
the Cu-0 bond length (0.05 A) and Cu charge (0.1e~) 
going from CU2O2 (Cu2+) to CU2O3 (Cu3+). 

We can attribute the relatively small charge difference 
between bulk and surface Cu ions (0.15 e~), to the afore- 
mentioned completion between Cu^+-like and Cu^^-like 
environments. To support further this interpretation we 
have calculated the charge difference between 2+ and 3+ 
Mn ions in MnO and Mn2 03 using the optimized data 
provided in Ref. [ij. Indeed, the resulting value of 0.38 
e~, by far smaller than the expected difference of 1 elec- 
tron which would result from a simplified ionic picture, is 
twice larger than the corresponding Cu-difference in our 
frustrated CuO/SrTiOs superstructure. 

Finalizing the discussion of the electronic structure, 
the gap of the ideal TET2 phase of 2.7 eV is strongly 
reduced when CuO is grown on SrTiOs: for a coverage of 
3.5 unit cells the gap formed by CuO states is reduced to 
about 0.6 eV, which for the total system (now including 
substrate states) is even more reduced to about 0.5 eV. 
The average local magnetic moment for Cu is about 0.75 
fiB^ which is enhanced in comparison to 0.63 fiB of the 
bulk TET2 phase. It should be noted that the present 
calculation for the surface system (CuO)[n]/SrTiO3(100) 
was done for ferromagnetic ordering (in order to make 
the HSE calculation feasible) whereas for the bulk TET2 
calculation the most stable phase was found for some 
particular antiferromagnetic ordering^. However, in Ref. 
l2| studying several magnetic orderings it was found, that 
the local Cu moment is rather insensitive to the specific 
alignment of spins. 



C. Summary 

Summarizing, our first principles study on the 
(CuO) [n] /SrTiOs (100) adsorbate system based on a hy- 
brid density functional theory approach describes and 
explains the recent experiments^ rather well, which cor- 
roborates our findings and predictions for the ideal bulk 
(without SrTiOa substrate) TET2 phased The analy- 
sis of our calculated results for structural and electronic 
properties indicates that the physical properties of the 
Cu-0 bonds are rather peculiar. In particular, we find 
that the enormous structural anisotropy of the Cu-0 sub- 
lattice -which determines the experimentally observed 
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FIG. 5: (Color online) Layer by layer total and l- 

decomposed charge of Cu atoms (sphere radius of 1.16 A) 
of (CuO)[3.5]/SrTiO3(100) in comparison to the correspond- 
ing values for bulk monoclinic CuO (tenorite structure) and 
the ideal bulk tetragonal TET2 phase. ^ The placements of 
the TiOa substrate and the vacuum regions are indicated. 



tetragonal syranietry- can be understood in terms of a 
layer dependent evolution of the Cu ionicity which in- 
creases progressively towards the surface. Ultimately, 
the local structural and electronic Cu-O environment ap- 
pears very frustrated as a results of the coexistence be- 
tween two concurrent states attributable to an in-plane 
Cu^+-like and out-of-plane Cu^+-like arrangements. 
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